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We study, in d-dimensions, the random walker with geometrically shrinking step sizes at each 
hop. We emphasize the integrated quantities such as expectation values, cumulants and moments 
rather than a direct study of the probability distribution. We develop a 1/d expansion technique 
and study various correlations of the first step to the position as time goes to infinity. We also show 
and substantiate with a study of the cumulants that to order 1/d the system admits a continuum 
counterpart equation which can be obtained with a generalization of the ordinary technique to obtain 
the continuum limit. We also advocate that this continuum counterpart equation, which is nothing 
but the ordinary diffusion equation with a diffusion constant decaying exponentially in continuous 
time, captures all the qualitative aspects of the discrete system and is often a good starting point 
for quantitative approximations. 
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I. INTRODUCTION 

An interesting extension of the ordinary random walk 
problem is a model in which the walker takes steps with 
decreasing size at each hop. There has been some interest 
among the physics community in the study of a particular 
system of this type during the past years namely the case 
for which the step lengths are given by A" , where n is the 
time step and A < 1. The system and its generalities are 
exposed in Q and for a one dimensional walker and 
its first passage properties are studied in Q , again for the 
one dimensional case. In this work we extend the study 
to d-dimensions and try to lay out important aspects 
within the room available. 

Unfortunately there is no immediate exact solution to 
the probability distribution for such a walker and conse- 
quently one is almost forced to invoke numerical simula- 
tions. To that end we tried to shy away from listing many 
probability density plots for various dimensions and for 
various values of A. Instead we focused on quantities 
which do not necessarily require a direct knowledge of 
the probability density, such as moments, cumulants and 
other interesting expectation values. 

In implementing the choice for hopping directions we 
use random unit vectors instead of picking orthogonal 
directions. This approach, as discussed in the text, keeps 
rotational symmetry for the probability distribution at 
all times. Using orthogonal directions for hopping vectors 
could also be interesting in of itself for other purposes but 
we do not pursue them here. 

One of the important aspects of the system at hand is 
a generalization of the shift symmetries of the ordinary 
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random walker; namely if one shifts time with one step 
one also has to shift the position variable by one step and 
sum among all possible directions. This symmetry is al- 
tered due to the geometric change in the step lengths, 
the translation of the position variable is now accompa- 
nied with a scale change of distances. We use this prop- 
erty to show that a recursive technique is available for an 
exact computation of the moments of the distribution, 
even when time becomes infinity where such an approach 
is not available for the ordinary random walker. This 
scaling property married with a, 1/d expansion approach 
proves to be very useful in computing various other quan- 
tities such as fractional power moments. 

We use the l/d expansion technique to study lower 
order cumulants of the probability distribution and con- 
trast them to those of a gaussian. We find that to order 
1 /d the probability distribution has vanishing cumulants 
except for the variance which is exactly what one finds 
for a gaussian. However the properly normalized vari- 
ances do not agree and this shows that the distribution 
is not exactly gaussian but very close to it. To that end 
we show that the 1/d approach allows for a continuum 
counterpart system -not a naive continuum limit- and the 
mentioned discrepancy is tied to the fact that for A < 1 
there is no continuum limit. Nevertheless we show that 
the continuum counterpart system is one for which the 
diffusion parameter decays exponentially with time. 

Since A < 1, the walker eventually comes to a halt in 
the infinite time limit and this means that there will be 
a non-zero probability density defined in a finite support 
characterized by the maximum distance the walker can 
take. This is in contrast to the ordinary random walker 
(A = 1) where the density vanishes in this limit. Due 
to this effect the choice for the first step (of length 1) 
is strongly correlated to the final position of the walker. 
We present the study of two such correlations. The first 
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one is the expectation value of the cosine of the angle 
between the first step and the position after infinitely 
many steps are taken. We present the numerical sim- 
ulations and compare it to the 1/d expansion method 
developed in the text. The second long term memory ef- 
fect we study is the probability, in the infinite time limit, 
for the walker to be along the -fz direction, given that 
after the first step the walker is at +z. This observable 
does not have al/ d expansion -it depends on ^fd-, how- 
ever using the continuum counterpart model as a starting 
point we find a phcnomcnological function which fits the 
numerical data rather well for all d and A. 



II. THE MODEL 

Consider a random walker in d-dimensions for which 
the step sizes decrease exponentially with each hop. Then 
if the walker starts at the origin, the random variable for 
its position is simply given by 



With the normalization factor given by 



N 



rN 



(1) 



Here f„ are random unit vectors obeying the following 
statistical rules 



(2a) 
(2b) 



Consequently we have (in any dimensions) the following. 



= 



(r^ ) = 



1-X 



(3a) 
(3b) 

(3c) 



An interesting property of the walker is that for A < 
1/2 the walker will fall out of range of the origin once 
it takes the first step. In this case the walker will be 
confined in a circular band with limits 



1-2A 
1- A 



< r < 



1 - A 



(4) 



The distribution function for rjv can also be expressed 
easily if we make use of the conservation of probability: 
if the walker is somewhere at step N -\- 1 then it must be 
at a neighboring point at the previous step. We therefor 
get 



P(f,iV+f) = -^ Y,Pir- ^^f-N.N) 



(5) 



(6) 



From our notation it is immediate that the probabil- 
ities to pick the directions are evenly distributed. This 
equation is nothing but a total sum over the possible 
neighbors normalized by the total number of neighbors. 
If one uses random unit vectors to pick hopping direc- 
tions, Af simply stands for the integral over the angular 
measure in a d-dimensional space. 

From Q it is clear that as ^ oo the probability dis- 
tribution approaches a time independent function which 
we will denote as Px{r). 



A. Rotational Invariance 

As far as rotational invariance is concerned, there is an 
important property of the walker we have presented. To 
expose this let us consider the model in 2-dimensions: if 
one uses orthogonal directions to pick the equal weight 
random hopping vectors, the distribution function Pxi^) 
will not necessarily have rotational invariance. This effect 
may be understood clearer if one considers A = where 
it is most pronounced. Since for A = the walker stops 
after the first step the probability distribution is, 

^[S{r^-x) + 5{f+x) + S{f-y) + 5{f+y)] . 

This has no continuous rotational invariance, only a dis- 
crete one. This effect should remain as we increase A. 
In fact one can hope to gain rotational invariance in the 
infinite time limit only as A — > 1. The slowing down 
of the walker breaks the rotational invariance if one in- 
sists on picking the random vectors along the orthogonal 
directions. 

As described in this peculiarity is related to the 
fact that the support of the distribution function is not 
a regular lattice: it is a collection of points that evolves 
in time in a rather complicated way. In studying the or- 
dinary random walker introducing a lattice makes sense 
since one knows that at each time step the walker resides 
on the backbone lattice and that in the infinite time limit 
rotational invariance is restored. However in the model 
presented here choosing lattice schemes for the choice of 
random hopping vectors will not recover rotational in- 
variance in full as time goes to infinity. Using random 
unit vectors to pick a hopping direction has a priori ro- 
tational symmetry. This is the reason why we have im- 
plemented the model with random unit vectors in the 
beginning of this section. On the other hand depending 
on the particular natural system where this model can be 
applicable, choosing a lattice scheme can be justified and 
could be of interest in of itself. In this work we will be 
working with the continuous random unit vector scheme. 
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III. GENERAL PROPERTIES OF Pa 

In this section we try to expose general properties of 
the probability distribution. In subsection A we study 
the Fourrier transform technique and present the solution 
for arbitrary d. In the next part we discuss the scaling 
property of the distribution and show that it allows for 
a recursive method to compute moments of the distribu- 
tion. Subsection C is reserved for the development of a 
1 /d expansion technique which we use in part D to study 
the variance, skewness and kurtosis of the distribution 
and contrast it to a gaussian. We also show in the last 
subsection that a continuum counterpart system is avail- 
able as the first term in a 1/d expansion of the master 
equation (O . The continuum counterpart system is not 
a naive continuum limit and has a diffusion constant that 
decays exponentially with time. 



A. Fourrier Space Solution 

One can try to make an explicit use of (|SJ) by intro- 
ducing the Fourier transform P\ of P\ which which will 
satisfy the following equation, 



Px{k,N+l)=M Px{k,N) J2 exp 



k ■ fjq 



(7) 



In view of rotational invariance we can apply the conven- 
tion of picking k along the z direction in a d-dimcnsional 
space. This will result in the following, 



P(fc, N + l) = P(k. N) f exp \-tX^kx] . (8) 

J LUd ^ 

Here x is the cosine of the angle between and the z 
direction. The measure dujd{x) is 



gid,k) 



1 



' r 

diVdix) 



exp(iA: • f) 



cos(fcx) 



(11) 



Where we again have made use of the convention of pick- 
ing k along z direction and in passing to the cosine func- 
tion we have used the evenness of P\{k). This results in 
the following list. 



(12a) 
(12b) 
(12c) 
(12d) 
(12e) 
(12f) 
(12g) 



Here Jj and ji are i'th Bessel and spherical Bessel func- 
tions respectively. This table amounts to the following 
expressions for general d > 1 



g{i,k) 


— cos(fc) , 


Q{'2,k) 


= Jo(k) , 


Q{?>,k) 


= jo{k) , 


Q{^k) 


= 2Ji(k)/k, 


g{5,k) 


= 3ji(fc)/fc , 


Q{Q,k) 


= 8J2(k)/k2 


g{7,k) 


= 15j2(fc)/fc2 



g{2n,k) = (2n- 2)!! 
g{2n+l,k) = (2n- 1)!! 



Jn-i(k) 

jn-l{k) 



(13) 
(14) 



With these constructs one has access to the usual ma- 
chinery of the Fourrier technique. For example the mo- 
ments of the distribution function in the position space 
will be given by coefficients of P(fc, N) in an expansion in 
k. There is however a more useful albeit indirect method 
which has its roots in the scaling property of the distri- 
bution. 



dijJd{x) = gd{x)dx = (l — z^) ^ dx , (9) 

the region of integration for x is [—1,1], and LJd is the in- 
tegrated measure. The formalism is valid for all d though 
one has to exercise attention for d = 2 and d ~ 1. For 
d = 2 we have x = cos{9) and du!2 = dO with 6 e {0, 27r}, 
the integrated measure is 27r. For d = 1 there is no in- 
tegral, only a sum over the discrete values of a; = ±1, so 
the normalization factor becomes 2 

The solution to ||SJ) is easily obtained remembering that 
if the walker starts at the origin we have Px{k,Q) = 1. 
Then by iteration we have the following, 

N 

Px{k,N)^Y[ g{d,\^k) , (10) 

with 



B. Scaling Property of P\{r) and Px{k) 

From H10|) it is clear that the distribution in Fourrier 
space enjoys the following scaling rule along with time 
translation 

Px{k,N + l) = Px{Xk,N)g{d,k) , (15) 
or in the infinite time limit 

Px{k) = Px{Xk)g{d,k) . (16) 

This suggests connections to various interesting phys- 
ical systems. For example taking the logarithm of H16|) 
we find the following 

V'(fc) = /i(fc) -f V(Afc) (17) 
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with ijj = In Pa (fc) and h = —\nQ{d,k). 

In form (|17l) is exactly like the free energy per bond 
acquired by an exact decimation of a thermodynamical 
system on a lattice with hierarchy . If the system ad- 
mits one relevant critical parameter tc, then for t = t — tc 
close to zero the free energy is (see for example and 
i) 

f{T)^g{T) + b-^ f{Xr) . (18) 

with b being the lattice scaling factor, V the effective 
dimension. Shrinking paramater A in this context is the 
slope of the linearized renormalization group transforma- 
tion. So the logarithm of the Fourrier transform of the 
random walker we study is like a thermodynamical sys- 
tem near criticality with zero effective dimensions, the 
critical point in this analogy being A: = for our sys- 
tem. This is somewhat expected since as k —^ we probe 
longer distance behavior of the distribution, in this region 
since the walker eventually stops at r^"-"-^ = (1 — A)^^ the 
distribution has to vanish for larger values. This hints at 
a slope discontinuity an hence at a possible critical point 
at r""^. 

It has also been shown in this context that there are 
connections to q-oscillators 8] and q-differcnces @. 
These analogies are interesting in their own rights. But 
in this paper we aim at a different direction and will not 
pursue them further. 

This scaling behavior of the distribution depends on 
the fact that the steps of the walker shrink geometri- 
cally with parameter A. This is better understood if we 
remember that position random variable obeys the fol- 
lowing, 

nv+i =fi + XpN ■ (19) 

That is if r/v+i is a random variable enough to describe 
the system at time A^-|-l so is pN to describe the system at 
time N . Their statistical independence is characterized 
by the discrepancy introduced by fi. This behavior of 
the position random variables should also manifest itself 
at the distribution level and we should have 

n(F,«+l)^(^)i:ft(l^,iv). (20, 

Obviously if A = 1 (|20|l and |SJl are identical, for A 7^ 1 
they also contain the same information but this is some- 
what clouded by the scale transformation. 

In this paper we will be studying infinite time limit 
properties of the distribution. In this limit the scaling 
property reads as 

^00 = ^1 -I- A Poo • (21) 

Here Poo is a random variable equivalent to Too to de- 
scribe the system, only it is a statistically different one 



where the difference is manifested in the random unit 
vector fi. This equivalency means that at the level of 
the probability distributions 

Which follows from (|23as N ^ 00. The form in 
have been presented in [2| for the one dimensional case. 

C. Scaling Invariance of Expectation Values 

The equivalency represented by H21|) also means exact 
equality if one considers expectation values, that is one 
must have for any function F 

{F{foo)) ^ {F{poo)) ■ (23) 

We can exploit this scaling behavior as follows. Let 
us identify the first step ri as the z direction in a d- 
dimensional space. In view of 1211 we have the following 
for any function F , 

{F{r)) ^Njdp^ Px{p) F{Xp + z) , (24) 

which is nothing but a direct consequence of (|22|l . From 
know on we are omitting the subscript 00 as it is clear we 
are working with the distribution function at the infinite 
time limit. Also at this point we would like to choose our 
normalization prescription. It is better for our purposes 
to work with a distribution normalized along the radial 
direction only 10] . With this convention wc have the 
following 

(F(f)) ^ f^dp p^-^ P,{p) F{Xp + z) (25) 

J UJd 

This final expression could prove more useful than a 
direct application of the Fourrier transform technique. 
One immediate consequence of (I25II is that it provides 
a recursive algorithm for computing the moments of the 
distribution function. To expose this let us start with the 
following, 

= f^dp P,{p) (1 + X'p' + 2Xpx) . 

J 

Since x has zero expectation value under the measure 
dijJd{x)/ujd we have 

{r')= J dpp'-'P,{p)il + X'p') , 
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but by our normalization of Px this equation is equivalent 
to the following, 

which in turn yields the correct result, 

If one chooses F = we get the following equation 



(1- A4)(r4) = 1 + 2A2 (r2) 



2 



where we have used the fact that has expectation value 
1/d under dijJd{x)/uJd giving in the correct result [TH 



1 



1 

TTa2 d 



(26) 



Higher moments get more complicated but are still 
computable with not so much effort. Therefor we have 
a useful recursive algorithm for finding all the even mo- 
ments of the distribution function. What we need to ap- 
ply the recursion are the general moments of the measure 
du)d{x), 



{x^ 
{x^-} 



dujd{x) 
a 

T{d/2)V{n- 



1/2) 



T{l/2)T{n + d/2) 
. 



(27a) 

(27b) 
(27c) 



D. 1/d Expansion of Expectation Values 



with {x"} = 1. For d ^ oo one has the following solution 



(2n- 1)!! 



(30) 



Of course in actual computations one must use the exact 
recursion 129|) . The point we wanted to expose is that 
this will allow for a l/d expansion of (|28l) in terms of the 
even derivatives of F{^J\ + -t- 2Xpx) with respect 
to X evaluated at a; = 0. The meaning of this is clear 
if we remember the meaning of the measure dujd{x) in 
@ : as d oo for all x except for x = this measure 
represents a number smaller than unity raised to a very 
large power. At a; = the measure is 1. Therefore as 
d ^ oo the normalized measure dujd{x)/uJd represents a 
properly normalized S function. 



lim — iLJ. —j. S{x)dx 

d^oo LUd 



(31) 



We therefor come to the following in which we use a 
condensed notation 



(32) 



ni=0 

where we have 

7.2n 



F^"nsi(r,0,A)) = 



(i2»iF(si(r,a;, A)) 



|x=o , (33a) 

Si(r,a;,A) = y^l + AV^ + 2\rx . (33b) 



The expression in (|32|l is not a complete expansion in 
1/d since one stih has to evaluate < F^"i (Vl + A^r^) >. 
But we can repeat the whole procedure between equa- 
tions H28() and (|32|l many more times. And the resulting 
form after k iterations is 



For more general functions F one could proceed as 
follows. We have the following direct from our formalism 



oo oo 



(Fir)) 



duJd 



dpp'^-^Px{p)F{^l + A2p2 + 2Xpx). 



(28) 

The complexity of this approach depends on the func- 
tion F. For general dimensions there is no great simpli- 
fication and one has to directly implement the formal- 
ism. Odd dimensions are better in view of the fact that 
the measure will not have a square root, particularly for 
d = 3 one has du)d{x) = dx, which tends to simplify the 
approach somewhat. On the other hand after a closer 
look at the moments of the measure du!d{x) in H27(l one 
sees the possibility of a 1/d expansion since one has 



(29) 



X (F2"---2»'=(sfc(r,0,A))) . (34) 
With the following terms 



5(ni,--- ,nk) = n o I 
-■- 2n, ! 



Sfe(r, x. A) 



1-A2fe 
1 - A2 



+ A^'^r^ -I- 2X^^-^rx 



^2n,j^2ni,...,2n,_i i^Sk{r,X,X)) 



|x=0 



(35) 



(36) 



(37) 
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Considering an infinite iteration of tlie procedure we 
arrive at the completed 1/d expansion 



ni=0 noo=0 



X F 



2ni,-" ,2n, 




1 



1-A2 



• (38) 



It is instructive to note tliat in this expression we eval- 
uate i^2ni,- - ,2noo a,t the value r = l/Vl — implying 
that they are evaluated bX d = co. The reason for this 
is the same as the large N argument for 0{N) field the- 
ories: if = oo only planar diagrams will contribute 
to the generator of connected green functions and all the 
higher correlations factorize in terms of the two-point ex- 
pectations. The correspondence to our model is that at 
d = oo, our d-dimensional random walker in the infinite 
time limit, is like a zero-dimensional 0{d) field theory. 
So when d = oo the only non-trivial moment is < > 
and all the others are given as 



lim (r^ 



So consequently we have 



lim (F(r)) = F 



A2 



p/2 



(39) 



(40) 



As it is implicitly assumed, for this expansion to work 
the F should be a smooth function. 



negative kurtosis excess means that the peak of the distri- 
bution is fiatter (platykurtic, suhgaussian) than the gaus- 
sian distribution for which C = 0. For A ^ 1 we recover 
a vanishing kurtosis excess as one would expect, since 
for this value of A the distribution is a gaussian in the 
infinite time limit. 



2. 1/d expansion of cumulants in d-dimensions 



In rf-dimensions, in our normalization scheme, one 

should actually pay attention to r"^^^ P\{r). This func- 
tion generally has non- vanishing odd moments {^r^"^^^. 
To obtain the kurtosis one needs n = 1 and n = 3 along 
with our exact results for ^r^)) and (r^)- Below, we 
present the results of a 1/d analysis up to and includ- 
ing order l/d^. The cumulants are defined as, 



C3 = 2(r)'-3(r)(r2) + (r3) , 

C4 = -6(r)Vl2(r2)(r)'-3(r2>' , 

The final form of cumulants are as follows 



A2 



C2 = 



dil + X") 



C3 = - 



A^(-7 + .5A-^ + .3A^) 
2d(l + A2)(l + A2 + A^) 

1 3A4(1 - A^ - A^) 

d? (1 + A2)2(1 + A2 + A4) ' 



(42a) 
(42b) 

(42c) 



(43a) 
(43b) 



E. Analysis of Lower Cumulants 

We do not know a priori the form of the distribution 
P\ and an exact treatment is not immediately evident. 
However, one can obtain valuable information for the 
form of the distribution by studying lower moments or 
cumulants. Two of the most important quantities are 
the skewncss and the kurtosis of a distribution which we 
will study below. 



1. Lower cumulants in 1 dimensions 

In one dimensions we have the following from previous 
sections 



C4 = 0+O(^), (43c) 
where we have normalized the cumulants as follows 



(44) 



2\"/2 ■ 



We would like to compare these cumulants to that of 
a gaussian in d-dimensions 



G(r) = Nd r"-^ exp 



2a 



(45) 



Here Nd is the normalization factor and cr^ = a^/d so 
that (r^) — 0-2 regardless of dimensions. Performing the 
same analysis for G we find the following, 



S = , 



A2 



2\2 



1-h A2 



(41a) 

(41b) 



Skewness S vanishes as one would expect for a symmet- 
ric distribution. However the kurtosis excess C < 0. A 



C3 



C4 



1 

2d 
1 
4d2 

+ 



(46a) 

(46b) 

(46c) 
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which are expansions of the results one can obtain in 
closed form for G. 

A comparison of cumulants in 143|l and H46|l reveals an 
exact match for A = 1 as one should expect. On the other 
hand, it is also evident that to order 1/d, the cumulants 
of Px have the form of the gaussian, albeit with a scale 
mismatch for the variance. That is, to order 1/d the 
distribution of the random walker we study is very close 
to a gaussian. This can be substantiated if we go back to 
the master equation |(SJ for Px and expand up to order 
1/d 



\2Arr 2 

PA(r,iV+l) = P;,(f,iV) + — — V2p;,(f) + --- , (47) 



2d 



where we have used 



d 



and introduced a length scale Sx via r — > Sxf. One can 
also introduce a time scale St so that one will end up with 
the following form 



dP 
'dt 



1 5x^ 
2d ~St 



(48) 



One can now try to do the usual practice of going to 
the continuum limit by taking Jcc — > along with (5i — > 
such that the ratio 6x^ /5t is the diffusion constant 



lim lim — — = D 

5x^0 &t^0 Ot 



(49) 



However there is a nuance; the factor A^^ with N = 
t/ 5t spoils this limit. So one must require the following 
as well, 



hm lim A^*/''* 

5*^0 A^l 



f{t/r). 



(50) 



Where the function / depends on how the limit is taken. 
The most natural choice for / that has all the flavor of our 
model is the exponential function fit/r) — exp(— i/r). 
That is we take the limit St ^ along with A ^ 1 such 
that 



'ln(A) 



(51) 



The complications that arise in this would be con- 
tinuum limit are related to the fact that in truth the 
random walker with shrinking steps does not walk on a 
static background of lattice points. The support of the 
distribution for any discrete time TV is a set that evolves 
in a complicated way everytime the walker takes a step. 



which prevents a continuum limit for A 7^ 1. We there- 
for arrive at the following 1/d continuum counterpart 

system 



dP 
'dt 



D 
2d 



exp(-t/T) V^P 



(52) 



This equation can better be written in terms of a dif- 
ferent time variable 



r = T[l-CXp(-i/T) 



(53) 



the limits t E (0,oo) corresponding to T e (0,1). And 
the equation becomes the common diffusion equation 



^ = ^V^P 
dT 2d 



(54) 



In 0, where the first passage characteristics of the 
random walker with shrinking steps were studied, is was 
shown that the continuum counterpart (|54|l captured all 
the qualitative aspects of the system with the equivalence 
of parameters r ^ — l/ln(A). 

We therefor have established that to 0{l/d) the dis- 
tribution function P\ is like a gaussian with a non- 
trivial rescaling of the variance, which is due to the non- 
triviality of the limit that resulted in the continuum coun- 
terpart equation. 

There are deviations in the cumulants in H43|) and 14t)|) 
starting at 0{l/d'^). The variance C2 > for all A and 
there is a critical value of Xvar which can make the 1/d^ 
term to vanish (but C3 survives it) 



v/33/5- 1 



The skewness given by 63 changes sign when the value 
of A becomes the square root of golden ratio. 



Xs 



Which means that to 0{l/d) and for A = Xskew the dis- 
tribution P\{r)r'^^^ is symmetric around < r >. 

Higher order expansions are in principle possible and 
could give further insight for Px but we do not pursue 
them here. 



IV. LONG TERM CORRELATIONS 

In this section we study some of the long term memory 
effects of the walker at hand. Since the walker ultimately 
comes to a halt due to the shrinking of the step sizes, it 
can not be very far away from the position of the first 
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step of unit length, which we can always choose to be 
the +z direction in a d-dimensional space. Thus the final 
position will have some memory of this first step unlike a 
normal random walker where the memory is washed out 
in the infinite time limit. 



A. Simple Examples: (z) and (2^) 

It is fairly straightforward to study quantities related 
to lower moments of the distribution. For example what 
would be the average position as time goes to infinity 
along the z direction given that the first step was along 
the positive z direction? Alternatively one could ask the 
same for the mean square displacement along the z di- 
rection. Applying the scaling property we get 



(^00) = 1 , (55) 

These are fairly easy to understand. The mean (zoo) 
would be 1 even for the ordinary random walk. But the 
fact that our walker eventually stops changes the behav- 
ior of (z^) which is expected to diverge as time goes to 
infinity, in the case of the ordinary random walk. The 
fact that (z^) ^ 1 as d ^ 00 is also expected: as the 
number of dimensionality grows the possibility to walk 
along a single direction decreases. In fact as d — > 00 all 
the subsequent walk following the first step can, to a very 
good accuracy, be considered to be happening along the 
directions orthogonal to z. 



B. The correlation {cos(S)) 

Our walker will not be infinitely far away from its first 
step. This, in turn, will introduce an angle correlation 
between the first step which is at z and the last step 
foo ■ An interesting quantity to study is the cosine of this 
angle, 



^ (iIm) ^ 

For the ordinary random walker this correlation should 
vanish in the infinite time limit since the walker would 
have spread all over the space equally likely. In this sub- 
section we will analyze this angle distribution for one and 
three dimensional cases and finally present a, 1/d expan- 
sion. 



1. One dimensional walker 

In one dimensions )12J the correlation we are after sim- 
plifies to 

{cos{0)) = (Sign(z)) (58) 

That is, we would be computing, in the infinite time limit, 
the probability that the particle is at z > minus the 
probability that it is at z < given that the first step is 
to the right. This will be related to the ultimate survival 
(remaining at z > 0) probability of a walker if the first 
step was to the right. 

It is not difficult to arrive at the following result after 
a direct implementation of the scaling property 

(Sign(z)) = / Px{z)dz (59) 
Jo 

At this point we would like to remind an important 
point about the limits of integration in general. Since 
the maximum distance the walker can travel is 1/(1 — A) 
the distribution will vanish for greater displacements . 
For A < 1/2 we have the following situation 

A < ^ =^ ^ > ^ (^'^^ 
which in turn means that 



A < - =^ (Sign(z)) = 1 (61) 

Thus for A < 1/2 the cosine of the angle between the first 
step and the last position are fully correlated, the reason 
is that the walker will fall out of range of the origin once 
it takes a step to the right. 

On the other hand as A —> 1 we have P\ and hence 
the correlation vanishes. The walker cannot remember it 
took the first step to the right as one should expect. 

The results of a numerical study of the problem are 
represented in Fig. We have also attempted an an- 
alytical study of for A = 2"^/™. It has been shown in 
01 and that for these values of A the distribution Pa 
are m'th order convolutions of a uniform distribution and 
hence a direct analytical evaluation is possible, albeit via 
symbolic programming packages. The squares in Fig. ^ 
denote results for to G A/". 



2. Three dimensional walker 

For a three dimensional walker we have to evaluate the 
following 

{cosm dxdu u^i^^ P^{u) (62) 

2 J |z-|-Au| 
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FIG. 1: Plot of ri{X)d for d = 1, defined in The back- FIG. 2: Plot of ri{X)d for d = 2. The dashed line represents 

ground points are results of a numerical simulation. The the 1/d expansion in the text. The other plot is the result of 

squares represent analytical results for A = 2"^'^™, with m a numerical simulation, 
from 1 to 30. 



The denominator can be expanded using the Legendre 
polynomials 



1 



Awl 



1=0 



1 



(63) 



with M> and m< are either 1 or Xu depending on which 
one is bigger. 



There is again a simplification if A < 1/2. In this case 
when u< = Xu the upper limit of integration on u is 1/A 
which is greater than the maximum range. So in this case 
any function of u under the integrand can be seen as an 
expectation value. There is a further simplification for 
d = 3. Since the angular measure is just dx we can use 
the orthogonality properties of the Legendre polynomials. 
This results in 



A < - for d = 3 
- 2 



3-4A^ 
3(1 -A2) 



(64) 



Unfortunately such a simplification is not present for 
A > 1/2. We present a computer simulation of (cos(6')) 
in Fig. 12 along with the result in H64|l . 
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FIG. 3: Plot of ri{X)d for d = 3. The dashed line represents 
the 1/d expansion in the text. The circles are representative 
points for the exact result H64|l obtained for A < 1/2. The 
other plot is the result of a numerical simulation. 



catcd rather fast and allows for an exact treatment only 
if A < 1/2. One can, on the other hand, use the 1/d 
expansion formalism we presented. Some algebra yields 
the following up to 1/ d^ 



3. d-dimensional walker 

In general dimensions one can still make use of the 
Legendre polynomial expansion, but this gets compli- 



(cos(0)) = V1-A2 



l + a(A)l + 6(A)l + 



(65) 

The first term is the d = oo value. As can be easily 
checked (cos(0)) =< z/r > and since in this limit corre- 
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Ultimate Fate to be in z > 
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FIG. 4: Plot of 'ri{\)d for d = 4. The dashed line represents 
the 1/d expansion in the text. The other plot is the result of 
a numerical simulation. 



lations factorize we can use < z >= 1 and r ^ = 



Vl — A^. The other terms are 



a(A) 



A^(1-2A^) 
2(1 + A2) 



(66a) 



Another interesting quantity is the probability p'^ , as 
time goes to infinity, to be at positive z values given that 
the walker took the first step along the +z direction. The 
reciprocal quantity p~ is to be at z < 0. These quantities 
are related by the conservation of probability, 



P 



P 

P" 



(Sign(z)) 



(69a) 
(69b) 



It is clear, due to the non-diiTerentiability of the signum 
function, that the \/d expansion procedure will not be 
directly applicable as in the previous section. Neverthe- 
less the behavior at d = oo is clear: the walker should 
remain at the positive z part of space. The large number 
of dimensions makes it impossible to choose many steps 
towards a particular direction. This is also immediate 
from the scale invariance of expectation values: we have 
(Sign(z)) = ({Sign(l -I- Arx)} > and d = oo means x = 
resulting in — — 1 and hence p~ = 0. We also 
expect, as A ^ 1, that p'^ = P^ = 1/2: the usual ran- 
dom walker spreads evenly in the infinite time limit no 
matter what the initial position is. Moreover, we again 
have a simplification for A < 1/2: if the walker takes the 
first step in the -t-z direction the origin and hence the 
negative z values are out of range. Thus, for A < 1/2 we 
have p~ = Q. 

The difficulty of a direct computation is revealed if we 
find the integral expression for , 



Q \4 

b{\) ^ — (16A8 + 12\^ - 13A4 - 29A2 + 15) 
8 



(66b) 



Higher order expansion is possible but we do not pursue 
it here. The series converges rapidly for this observable. 
We demonstrate this fact by plotting the following quan- 
tity both for a set of numerical simulation and the 1/d 
expansion 



77(A)d^(cos(0)),-(cos(0)),^ 



d—oo 



(67) 



In d = 2 the agreement is marginal and one has to 
proceed to a higher order expansion here. For d = 3 
the agreement is good and already from d = 4 we have 
considerable convergence. 

As one can see ry(A) changes sign meaning that at this 
particular value of lambda the observable (cos 9) is given 
by its d = cxD limit. These critical values can be found 
from our 1 / d expansion and they are given as follows 



P = 



dr 



l/A 



Pxir) 



^^^^"^ dujd{x) 



(70) 



Since the region of integration on r do not cover the full 
support of the distribution we can no longer apply our 
1/d program. Nevertheless the integral over x can be 
carried out exactly and this results in the following, 



2p' 



drr'^-^Px{r) 



l/A 



1-6 



d- 1 



(71) 

Here B{x^ a, b) is the regularized incomplete beta func- 
tion. 

We can gain further insight into the d behavior of p^ 
if we are content with the continuum counterpart model 
we presented in (|54|l for which the solution is simply the 
solution of an ordinary random walker with t replaced 
with T in (|53ll . In this scheme we can assume that the 
walker initially starts at f*o = ZoZ and find the probability 
to be at positive z values in the infinite time limit. A 
simple computation yields. 



1 



1 - 



56d 



(68) 



Pcont 



l + Erf 



/dz|_ 
2Dt 



(72) 
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other dimensions out. This will result in a one dimen- 
sional distribution which will be the same as the distri- 
bution of a one dimensional walker taking steps of length 
A"cos(6'„) at each time step n, with cos(6'„) distribution 
given by dujd{x). Every observable is directly related to 
the random walker with shrinking steps in one dimen- 
sions, for example one will have < >d=< >d=i /d. 
Thus most of the d and A dependence might come from 
(|72|l . Pursuing this observation and the parametric anal- 
ogy between r and A 



ln(A) 



we can follow a more phenomenological approach, 
can try to fit to the following estimate function 



We 



FIG. 5: Simulation results of for various dimensions. A 
range is upto 0.99. 




FIG. 6: Percent error between the simulation data and the 
fit function in 1731 1 with parameters given in TABLE I. 



Pestm 



lAid) 



1 + Erf 



m 

V2 



v/-rfln(2A- i; 



(73) 

which has all the global properties of p+ mentioned at 
the beginning of this section. 

In TABLE I we present the results of fits of p+ obtained 
from numerical simulations to the function in (|73|l . 
represents the chi-squared value of the fit and c is its 
correlation factor. As can be seen the fit is generally 
good throughout the range of d. As expected the fit has 
an increasing accuracy with increasing d. 

In Fig. |S1 one can find the simulation results for 
for various dimensionalities. In Fig |^we show the per- 
centage error ,as a function of A, between the simulation 
results and the fit function with parameters given in TA- 
BLE I. As can be seen the agreement is good for all values 
of A. 



d 


A(d) 


B(d) 


x' 


c 


1 


1.00305 


0.974315 


4.876 10"* 


0.999674 


2 


1.01298 


0.976132 


6.222 10"'' 


0.999708 


3 


1.00667 


1.003400 


6.491 10~* 


0.999596 


5 


1.00330 


1.015200 


3.378 10"'' 


0.999702 


7 


1.00183 


1.01426 


1.747 10~* 


0.999800 


9 


1.00107 


1.01129 


8.787 10"^ 


0.999875 



where Erf is the error function. Apparently the depen- 
dence on d is not what one could recover via a 1/d expan- 
sion. However, let us remember, from previous sections, 
that to order l/d, P\ is very close to its continuum coun- 
terpart system. This, hints at the possibility that the 
continuum counterpart result could be a very good ap- 
proximation for any d. This proposition can be expanded 
a bit better as follows. If one wants to compute p'^ one 
can just as well proceed first by computing the proba- 
bility distribution to be at (z, z + dz) by integrating the 



TABLE I: Fit parameters of the simulation data to the esti- 
mate function in II73|I . 



V. CONCLUSIONS AND FUTURE 
DIRECTIONS 

We have tried to expose some of the important prop- 
erties of a random walker with step sizes shrinking geo- 
metrically with each step. The most useful property in 
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general is the scaling property of the probability distri- 
bution when one increments the time step by one. This 
allowed for a \/d expansion approach the same way it 
would for the ordinary random walker. However, in the 
model studied this approach remains useful in the infinite 
time limit where the distribution goes to a time indepen- 
dent function. 

We have presented some interesting correlations of the 
position at large times to the first step. To that end we 
showed one correlation with a nice 1 /d expansion and one 
without. The continuum counterpart system is shown to 
capture all the qualitative aspects of the system and in 
general it proves to be a good starting point to approxi- 
mate the quantitative properties as well. 

As in most statistical systems one, two and three di- 
mensional systems are more interesting than higher di- 
mensions where the 1/d approach -actually reminiscent 
of a mean field approach- becomes more and more re- 
liable. One hint for this is already apparent in FIG 5 
and TABLE I, where the curves and fit parameters are 
characteristically different than the ones for higher di- 
mensionalities. The common topological differences be- 
tween two and higher dimensions remain. Consequently, 
for low dimensions, a more incisive approach to the sys- 
tem possibly including the direct study of the probability 



distribution is important and is left for future work. 

Another possibly interesting avenue to explore is to 
use lattice schemes to pick hopping directions. Similarly, 
this should matter only for low dimensions since the mean 
field should set in regardless for higher dimensions. 

We have also left out the first passage properties of 
the system at hand since it requires a study of the dis- 
tribution for finite times. This could also prove to be of 
interest as a topic for further research. The continuum 
counterpart system should prove to be a nice starting 
point, for example, for the survival probability in z > 0. 
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